  figure(3);
 count = 0; 
 h=8; 
 
plot(nn(floor(win/2)+h:floor(win/2)+size(MSFE_var{h}+h,2)+h-1),[ prctile(MSFE_ar{h},[16 50 84],1)'],'r'); datetick('x','yyyy','keepticks'); 
hold on
plot(nn(floor(win/2)+h:floor(win/2)+size(MSFE_var{h}+h,2)+h-1),[ prctile(MSFE_ar{h},[50 ],1)'],'k'); datetick('x','yyyy','keepticks'); 
 hold on
 set(gca, 'XTick', 695496:7305:732021);
 set(gca, 'XTickLabel', anni(:,1));  
 anni=datevec(695496:7305:732021);
plot([datenum(1914,9,15) datenum(1914,9,15)], [0 5000]); axis([695405 732021 0 5000] )
plot([datenum(1931,6,15) datenum(1931,6,15)], [0 5000]); axis([695405 732021 0 5000] )
plot([datenum(1941,9,15) datenum(1941,9,15)], [0 5000]); axis([695405 732021 0 5000] )
plot([datenum(1951,3,15) datenum(1951,3,15)], [0 5000]); axis([695405 732021 0 5000] )
plot([datenum(1971,9,15) datenum(1971,9,15)], [0 5000]); axis([695405 732021 0 5000] )
plot([datenum(1983,12,15) datenum(1983,12,15)], [0 5000]); axis([695405 732021 0 5000] )
hold off
title('FORECASTING INFLATION AT 8-QUARTERS HORIZON WITH OUTPUT GROWTH')
